Members of the Anopheles gambiae complex are the most important malaria vectors in the African continent. It is a group of seven morphologically undistinguishable species of tropical mosquitoes (Service, 1985; Hunt et al., 1998). Within this species complex the most important malaria vectors are: Anopheles gambiae (Giles, 1902) and An. arabiensis (Patton, 1905) which are distributed over 70% of the Sub-Saharan Africa with An. gambiae being restricted to the more humid and forested localities and An. arabiensis being distributed over the dry savanna and semi- arid parts of the African continent (Service, 1980; Bryan, 1983; Lindsay et al., 1998). The two species are more adapted to the human environment;they are sympatric and synchronic over most of their distribution range (Petrarca et al., 1998). Female An. gambiae complex mosquitoes have slender body, (5mm) in size with three sections: head, thorax and abdomen. The morphological description of the complex had been given by Evans (1938), cited in Gillies and De-Meillon (1968).
Due to the observed differences in vectorial capacities and behavioural habits of the siblings in species complexes, efforts have been geared towards identification of the siblings of mosquito species. One of these efforts is the use of multivariate methods such as morphometrics (Coluzzi, 1964; White, 1977; Lambert and Coetzee, 1982; Service, 1988; Schmidt et al., 2003). Morphometrics is the field concerned with variation and changes in the body form or size of an organism. It transforms the complex body forms into quantitative series of numbers that can be analyzed and used for comparisons between the different forms (Daly, 1985).
Since the discovery that An. gambiae is not a single species but a complex of sibling species, many researchers attempted to find morphological variations between the siblings of the complex. These attempts had been done to some extents by Ribbauds (1944), Muirhead-Thomoson (1951) and Coluzzi (1964). Subsequently, many researchers have attempted to find morphological differences between the three fresh water species: An. gambiae, An. arabiensis and An. quadriannulatus (Ismail and Hammoud, 1968; Zahar et al., 1970; White and Muniss 1972; Reid, 1975a) without success. A subspecies An. quadriannulatus davidsoni, from Cape Verde Island was described based on morphological criteria only by Ribeiro et al. (1979). Ribeiro (1980) proposed a phylogeny for the complex, An. gambiae using the morphological data from An. quadriannulatus davidsoni and published data for other members of the complex. Coetzee (1989) carried out a morphometric analysis on all life stages of An. gambiae, An. arabiensis, An. quadriannulatus and An. merus occurring in southern Africa. He revealed that the length of the hind leg pale band at the junction of tarsomeres 3 and 4 is a good character for grouping An. gambiae/ An. arabiensis and An. quadriannulatus/ An. merus and he was able to separate An. quadriannulatus/ An. merus in the base of palpal index. He found that separation of individuals An. gambiae from An. arabiensis was not reliable. Petrarca et al. (1998) carried out a morphometric multivariate analysis on field and laboratory specimens of An. gambiae and An. arabiensis from different areas of Sub-Saharan Africa and he found that all the measures were significantly larger for An. arabiensis.
Separation between the species of An. gambiae complex is much less reliable owing to the existence of a considerable overlap between the different character distributions (Ribeiro et al., 1979; Bryan 1980; Bushrod, 1981).
In the present study we carried out a morphometric analysis on field specimens of An. arabiensis and An. gambiae (previouslyidentified by molecular polymerase chain reaction (PCR) techniques) in an attempt to find reliable morphological features for distinguishing between the two species.
1 Materials and methods
1.1 Collection site
Adult female Anopheles mosquitoes were collected from Wau town, Republic of Southern Sudan which is situated in the eastern part of the African continent. Wau is located in mid western part of Bahr El Ghazal region between 7°: 42Ë N and 28°: 283Ë E. The region is an area of natural swamps and ironstone plateaus characterized by equatorial (humid-tropical) climate. The climate in this region is marked by high temperature, high rainfall and very high humidity.
1.2 Mosquito collection and preservation
In door resting wild adult Anopheles mosquitoes were collected by hand capture using sucking tube (aspirator) during the rainy season 2010. Anopheles mosquitoes collected were preserved individually in 70% ethanol and kept at -20℃ for subsequent processing in the laboratory. Samples were transported to the laboratory well protected to minimize any damage. The processing of the materials for this study was carried out at the Department of Zoology, Faculty of Science, University of Khartoum, Sudan.
1.3 Identification and mounting offemales An. gambiae complex
An. gambiae species complex were identified to the complex level using morphological identification keys described by Gillies and De-Mellion (1968) and Gillies and Coetzee (1987). Subsequently, females An. gambiae complex were dissected for morphological and molecular identification. For morphological analysis, permanent slides were made for detailed examination of external body structures of females An. arabiensis and An. gambiae. The head, wings and legs were carefully separated from the mosquito body and mounted in Puri's mounting medium on glass microscopic slides as described by WHO (1975) with minor modifications. The mosquitoes were mounted directly into Puri's mounting medium without prior clearing because gum-chloral mountant continues to clear specimens after mounting. The corresponding carcasses were preserved in 70% alcohol, kept at -20℃ and used for molecular identification.
1.4 Selection of morphological characters for morphometrics analysis
For terminology and morphological characters and abbreviations of palps, antennae, legs and wings, the nomenclature of Evans (1938) and Gillies and De-Mellion (1968) were adopted, except for wing spots, the nomenclature adopted by Gillies and Coetzee (1987) was followed.
A preliminary list was prepared from all morphological characters of the adult female An. gambiae complex include characters which were used by the previous workers in an attempt to discriminate between the members of the complex. Out of this list, all characters of external body partsthat fixation and transport could not have changed were selected. Sixty nine characters (2 non metric and 67 metric) were finally retained and were measured on 53 females An. arabiensis and 35 females An. gambiae. The measured characters include 2 characters associated with antennae, 11 characters with palps (Table 1), 14 characters with wings (Table 2) and 14 characters of each fore, mid and hind legs (Table 3).
Table 1 Mean, standard error and sample size of 2 antenna and 11 palp morphometric characters of females An. arabiensis and An. gambiae
|
Table 2 Mean, standard error and sample size of 14 morphometric wing characters of females An. arabiensis and An. gambiae
|
Table 3 Mean, standard error and sample size of 14 morphometric fore leg, mid leg and hind leg characters of females An. arabiensis and An. gambiae
|
1.5 Morphological characters measurements
Measurements were done by the projection method (Zahar et al., 1970). Selected morphological characters of females An. arabiensis and An. gambiae were measured using a Wild MII binocular calibrated compound microscope fitted with a 1.25× phototube. Characters were measured using 3.5× and 40× objective lens depending on the size of the measured characters. The slide mounted specimens were projected with a phototube on a microscopic field using 8× eye piece. After an excellent view of the projected image of the measured character was obtained, the projected image was drawn and subsequently measured to the nearest half millimeter. Then the lengths were calibrated to the real lengths using a micrometer stage (1mm-Erma, Tokxc).
1.6 Statistical analyses
The soft-ware computer program package SPSS® version 16.0 for Windows was used for statistical aspects of the morphometric analysis.The measurements were not transformed to ratios so as to preserve the possible influence of differences in the body sizes of the two species. Box plots were used to check for the presence of outliers (Garson, 2012). Then the measured data were subjected to univariate and multivariate statistical analyses.
Univariate statistics like calculation of descriptives, Students t test (t-test) and Box and whisker plots were used.Descriptives analysis involvedcalculation of mean size measurements plus or minus their standard error and sample sizes. A t-test was carried out to test the significant differences of measured morphological characters among An. arabiensis and An. gambiae populations (Park, 2009). Box and whisker plot were used to summarize the univariate data.
Multivariate statistics such as principal component analysis and discriminant function analysis were used. Principal component analysiswas conducted to simplify subsequent analysis of the data. Thetest is essentially a method of data reduction that aims to produce a small number of derived variables that explain most of the variance and can be used in place of the larger number of original variables (Pimentel, 1992). The analysis generates a set of principal components by weighting all the available variables. The first component explains the most variation; the second explains the next most variation, and so on. The amount of variance captures by one component represents by Eigen value. Investigation of the first few components will show which variables contribute most to the variations between individuals (Dythan, 2003).
Discriminant function analysis concerns with classification and aims to obtain a small number of useful discriminating variables (Pimentel, 1992). Discriminant function analysis was carried outbased on the results of the principal component analysis and t-test to assess the degree of discrimination among adult female An. arabiensis and An. gambiae.
2 Results
2.1 Morphometrics investigations
Tables 1, 2, 3 summarize the mean, standard error and sample size of the measured antenna, palp, wing and legs characters on females An. arabiensis and An. gambiae.
2.2 Student t-test
The t-test showed significant variations (P<0.05) between the sympatric populations of An. arabiensis and An. gambiae in 17% (13/69) of all the measured morphological characters (Table 1, 2, 3). Most of the significantly different characters (9/13) were larger in An. arabiensis individuals. Variations in the measured characters show promise for a discriminant function analysis in an attempt to distinguish between the two fresh water species.
2.3 Principal component analysis
Principal component analysis was carried out to simplify the analysis of Tables (1, 2, 3).The results of this analysis as presented in (Table 4) showed that three components with Eigen values more than 1 explained the data. Component one contributed the highest Eigen value of 0.52 with a percentage variance of 83.2%. It was largely influenced by the measurements of the mosquito body size. Component 2 and 3 showed Eigen values and percentages variance of 0.04, 6.3% and 0.02, 3.2%, respectively. They were mostly influenced by the variables measures of the body shape. These variables were less important due to the lower percentages of variance contributed by their components. Together, component 1, 2 and 3 account for 92.7% of the total variance. All the variables positively load on the first component except the: palpal index, humeral pale spot, pale band at the joint of tarsomeres 1and 2 of the hind leg and pale band at the joint of tarsomeres 3 and 4 of the hind leg.
Table 4 Results of principal component analysis: contribution of 37 variables to the first three principal components calculated from 53 females An. arabiensis and 35females An. gambiae
|
The inspection of principal component 1 showed that individuals of An. arabiensis and An. gambiae were not significantly different (P>0.05). Box and whisker plot (Figure 1) summarize the scores of component 1.
Figure 1 Box and Whisker plots of the scores of principal component 1 (83.2% explained variance) of females An. arabiensis and An. gambiae populations. Value zero on the Y-axis is the grand centroid (overall mean of the components scores) (P>0.05)
|
2.4 Discriminant function analysis
Discriminant function analysis was carried out using two different sets of variables. First, the analysis was applied on the characters linked with the mosquito body size. These characters comprised the five most influential characters for component 1, derived by the principal component analysis (Table 4). The characters used were: hind leg length, mid leg length, fore leg length, wing length including the fringe and wing width including the fringe. When these characters were subjected to discriminant function analysis, one significant function was derived. A P-value (Wilks Lambda test) bigger than 0.0001, indicates that the group centroids of the 2 species were not significantly different. The canonical correlation coefficient of the function was 0.401 and the related Chi-square (χ2) =6.766, this indicates a low correlation between the discriminant function and the original variables. The Chi-square was not significant (P>0.0001, df=5), indicating populations with no definite differences between the two species. Table 5 shows the standardized, unstandardized and loading coefficients for the discriminant function. The variables that mostly influenced the separation of the two species were (with a decreasing order of magnitude): wing length including the fringe, fore leg length and mid leg length. However, this discrimination was not a complete one, since only 67% of the original grouped specimens and 58% of the cross validated ones were correctly classified (Table 6).
Table 5 Canonical discriminant coefficients and loadings for the discriminant function. The analysis applied on the characters linked with mosquito body size to discriminate between females An. arabiensis and An. gambiae
|
Table 6 Leave-one-out cross validation for all specimens used in the discriminant analysis of morphomretric measurements (characters linked with mosquito body size) of females An. arabiensis and An. gambiae
|
Second, the 13 characters that showed significant differences (using t-test) between the two populations of females An. arabiensis and An. gambiae (Tables 1, 2, 3) were subjected to the discriminant function analysis. One significant function was derived. A P-value (Wilks Lambda test) smaller than 0.0001, indicates that the group centroids of the 2 species were significantly different (P<0.0001). The canonical correlation coefficient of the function was 0.869 and the related Chi-square (χ2) =48.603, df =13, P =0.0000. The function explained the total variance among the two species, with (with a decreasing order of magnitude) sector pale spot, tarsomere 4 of the fore leg, pre- sector pale spot, tarsomere 5 of the hind leg and the tarsus of the fore leg having high discriminant loadings (Table 7).
Table 7 Canonical discriminant coefficients and loadings for the discriminant function. The analysis applied on the characters selected by t-test to discriminate females An. arabiensis and An. gambiae
|
The means of the discriminant scores (group centroids) of An. arabiensis and An. gambiae were 2.123 and -1.388, respectively. The difference between the centroids of the two species was highly significant (P<0.0001). However, this discrimination was a complete one, since only 85.2% of the original grouped specimens and 72.7% of the cross validated ones were correctly classified (Table 8), as indicated in Figures 2.
Table 8 Leave-one-out cross validation for all specimens used in the discriminant function analysis of morphometric measurements on females An. arabiensis and An. gambiae
|
Figure 2 Characters showed significant differences between females An. arabiensis and An. gambiae using t-test (P<0.05). These characters were obtained 85.2% discrimination when applied to discriminant function analysis
|
The two species were previously identified by the molecular polymerase chain reaction (PCR) techniques. Using this morphometric method, 11 (20.8%) of the 53 An. arabiensis were misclassified as An. gambiae. Of the 35 An. gambiae 2 (5.7%) were misclassified as An. arabiensis. The mean error rate was 13.3%.
Figure 3 shows the scatter plot of tarsomere 4 of the fore leg against the sector pale spot of the wing, the two variables that mostly influenced the separation of the two species. The scatter plot showed separation between most individuals of the two species.
Figure 3 Scatter plot of the length of tarsomere 4 of the fore leg against the length of sector pale spot of females An. arabiensis and An. gambiae. The two characters mostly influenced the separation of the 2 species as revealed by 85.2% discrimination
|
3 Discussion
Measurements of various morphological characters of females An. arabiensis and An. gambiae from Republic of Southern Sudan were in accordance with the published data of the same species from other sites along the distribution range of the two species (Coluzzi, 1964; Petrarca et al., 1998; Adeleke et al., 2008).
The results of the principal component analysis confirmed the report of (Dythan, 2003) that the individuals of females An. arabiensis and An. gambiae differed mainly in the mosquito body size measurements. The morphometric analysis (discriminant function analysis and t-test) showed that females An. arabiensis and An. gambiae were not significantly different in the body size measurements. In contrast Petrarca et al. (1998) found that An. arabiensis had mean body size greater than An. gambiae.
Using the significantly different characters between the two species, as revealed by t-test- the discriminant function analysis revealed differentiation between the two species to a confidence level approaching 85.2%. Characters that mostly influenced differentiation were: sector pale spot, tarsomere 4 of the fore leg, pre-sector pale spot, tarsomere 5 of the hind leg and tarsus of the fore leg. In other groups of insects a morphological discriminability bigger than 75-80% is sufficient to propose the separation to species of populations of differentiated morphs in dragon flies (Carrison, 1992).
The best morphological character selected by the morphometric analysis for discrimination of adult females An. arabiensis and An. gambiae was the sector pale spot of the wing. The spot vary significantly between the two species. Petrarca et al. (1998) stated that this spot is important for discriminating the laboratory strains of An. arabiensis and An. gambiae as observed by Coronel (1962) but at the same time the character failed to discriminate the field strains of the two species.
In this morphological investigation, all the individuals of An. gambiae and An. arabiensis examined for the number of pale bands on the palpus were 3-banded. The same finding has been reported by Gillies and De-Mellion (1968) and Bryan (1980). The palpal index was used by Coluzzi (1964) and Bryan (1980) in separating the salt-water species from fresh-water species of the An. gambiae complex. The results of the morphometrics analysis indicated that this character did not discriminate the two species, as observed by Petrarca et al. (1998).
Coetzee (1986, 1989) examined the use of hind leg banding patterns for identifying members of the An. gambiae group of mosquitoes. They showed that the pale band at the apex of hind tarsus 3 and the base of hind tarsus 4 was separated An. gambiae and An. arabiensis from An. merus and An. quadriannulatus. In our study, the character failed to distinguish the sympatric populations of An. gambiae and An. arabiensis.
Coluzzi (1964) stated that the number of coeloconic sensillae of the antenna might prove useful for the differentiation An. gambiae and An. arabiensis in East Africa. Also the partially discrimination value of the number of coeloconic sensillae is confirmed by the observations of other authors (Ismail and Hammoud, 1968; Petrarca et al., 1984; Coetzee, 1989). In our study, the mean number of antennal coeloconic sensillae did not vary significantly between the two species as pointed out by Petrarca et al. (1998).
Adeleke et al. (2008) carried out a morphometric analysis of An. gambiae complex in Abeokuta, Metropolis. The results suggest that the antennal length and wing length may be of significant value in separating An. gambiae and An. arabiensis. This results contrast sharply with the finding of our study that the two characters did not vary significantly between the two species.
The morphometric analysis conducted in this study confirmed the existence of variability within several characters among the populations of An. arabiensis and An. gambiae. This study suggests that a considerable degree of differentiation for some metric morphological characters exists among the females of the two species. More studies are needed on morphometric studies on An. gambiae complex so as to document as many discriminating characters as possible. Further studies are still recommended by comparing the measured morphological characters with genetic composition in order to give valuable clues for identifying the sibling species.
Authors' contributions
Hamza A.M. participated in field sampling, performed the technical work, analyzed the results and drafted the manuscript. Abukashawa S.M. contributed to the reading and revising the draft manuscript. El Rayah El A. supervised the research group and participated in revising the draft manuscript. All authors read and approved the final manuscript
Acknowledgements
The authors thank Dr. Noteila Mustafa Khalid, Department of Zoology, Faculty of Science, Khartoum College of Medical Sciences for her collaboration and assistance in the statistical analysis of the results. We thank Mr Abd El Gader Morwak, Ministry of Health, Western Bahr El Ghazal State, Republic of Southern Sudan for his assistance in the collection of mosquito samples from Republic of Southern Sudan. We thank Department of Zoology, Faculty of Science, University of Khartoum for the institutional support and use of the facilities of the Zoology Laboratory. This study was supported by a grant from Kassala University.
Adeleke M.A., Mafiana C.F., Idowu A.B., Amusan A.A., and Adekunle M.F., 2008, Morphometric studies on Anopheles gambiae complex (Diptera: Culicidae) in Abeokuta, southwest Nigeria, Acta Entomolegica Sinica, 51(12): 1289-1292
Bryan J.H., 1983, Anopheles gambiae and Anopheles melas at Brefet: The gambiae and their role in malaria transmission, Ann Trop Med Parasitol,77: 1-12
Bryan J.H.,1980, Use of palpal ratio and the number of pale bands on the palps in separating Anopheles gambiae Giles s. s. and Anopheles melas Theobald (Diptera: Culicidae), Mosquito Systematics, 12: 155-163
Bushrod F.M., 1981, The Anopheles gambiae Gilescomplex and Bancroftian Filariasis transmission in a Tanzanian coastal village, Ann Trop Med Parasitol, 75: 93-100
Carrison R.W., 1992, Using ordination methods with geographic information species resolution in a partially sympatric complex of Neotropical Tramena dragon flies (Odonata: Libellulidae), In. J.T. Sorensen and R. Foottit (eds.), Ordination in the study of morphology, evolution and systematic of insects, Elsevier, Amsterdam, The Netherlands, PP 223-340
Coetzee M., 1986, Practical use of the hind leg banding patterns for identifying members of the Anopheles gambiae group of mosquitoes, Mosquito Systematics,18: 134-138
Coetzee M., 1989, Comparative morphology and multivariate analysis for the discrimination of four members of the Anopheles gambiae group in southern Africa, Mosquito Systematics, 21 (2): 100-116
Coluzzi M., 1964, Morphological divergence in the Anopheles gambiae complex, Riv Malariol, 43: 197-232
Coronel L.T., 1962, Morphological variation in Anopheles gambiae Giles, WHO/ Mal/, WHO Offset Publications, Geneva, Switzerland, 328: 6-10
Dythan C., 2003, Choosing and Using Statistics: A Biologist/s Guide, 2 Malden, Massachusetts, USA: Blackwell Publishing
Evans A.M., 1938, Mosquitoes of the Ethiopian Region, British Museum, National History
Garson G.D., 2012, Testing statistical assumptions, Statistical Associates Publishing, Blue Book Series, Asheboro, USA, PP 30-31
Giles G.M., 1902, A handbook of the Gnats or Mosquitoes, 2nd edition. London: John Bale, Sons and Danielsson
Gillies M.T. and Coetzee M., 1987, A supplement to the Anophelinae of Africa South of the Sahara (Afrotropical region), South African Institute for Medical Research, Johannesburg, South Africa
Gillies M.T. and De-Mellion B., 1968, The Anophelinae of Africa South of the Sahara, South African Institute for Medical Research, Johannesburg, South Africa
Holstein M., 1952, Biologie d' Anopheles gambiae, Monogr, Ser WHO, 9: 176, (Rome), 13: 421-426
Ismail I.A.H. and Hammoud, E.I.H., 1968, The use of coeloconic sensillae on the female antenna in differentiating members of Anopheles gambiae Gilecomplex, Bull WHO, 38: 814-821
Lambert D.M. and Coetzee M., 1982, A dual genetic and taxonomic approach to the resolution of the mosquito taxon, Anopheles (Cellia) marshalli (Culicidae), Systematic Entomology, 7: 320-332
Lindsay S.W. and Mortens W.J., 1998, Medical Biology, Climate and Malaria, Bull WHO, 76: 33-45
Park H.M., 2009, Comparing Group Means: T-tests and One-way ANOVA Using STATA, SAS, R, and SPSS, Working Paper, The University Information Technology Services (UITS), Center for Statistical and Mathematical Computing, Indiana University, http://www.indiana. edu/~statmath/stat/all/ttest
Paterson H.E., 1964b, Salt-water Anopheles gambiae on Mauritius, Bull WHO, 31: 635
Patton W.S., 1905, The culicid fauna of the Aden hinterland, their haunts and habits, J. Bombay Nat Hist Sot, 16: 623-637
Petrarca V., Carrara G.C., Deco M.A., and Petrangeli G., 1984, Osservazioni citogenetiche e biometriche sui membri del complesso Anopheles gambiae in Mozombico, Parasitologia, 26: 247-259
Petrarca V., Sabatinelli G., Toure Y.T., and Di Deco M.A., 1998, Morphometric multivariate analysis of field samples of adult Anopheles arabiensis and Anopheles gambiae s. s (Diptera: Culicidae), J Med Entomol, 35(1): 16-25
Pimentel R.A.,1992, An introduction to ordination, principal components analysis and discriminant analysis, In JT Sorensen, R Foottit (eds), Ordination in the Study of Morphology, Evolution and Systematics of Insects, Elsevier Science Publisher, Amsterdam, PP 11-28
Reid J.A., 1975a, Pupal differences between species A and B of the Anopheles gambiae group from Kisumu, Africa, Mosquito Systematics, 7: 1-7
Ribbands C.R., 1944, Differences between Anopheles melas (An. gambiae var. melas) and Anopheles gambiae, Ann Trop Med Parasitol, 38: 85-99
Ribeiro H.H.C., Ramos C.A., and Capela R.A.,1979, Description and biometric study of Anopheles (Cellia) quadriannulatus davidsoni ssp, a seventh member of the Anopheles gambiae Giles complex (Diptera: Culicidae) endemic to the Cape Verde Archipelago, Garcia de Orta Ser Zool, 8: 75-88
Ribeiro H., 1980, A biometric study of the taxonomy of the Anopheles gambiae Gilescomplex (Diptera: Culicidae), Garcia de Orta Ser Zool, 9: 139-154
Schmidt E.R., Foley D.H., Burgoro H., and Bryan J.H., 2003, A morphological study of the Anopheles punctulatus group (Diptera: Culicidae) in the Salomon Islands, with a description of Anopheles (Cellio) irenicus Schmidt, spn, Bulletin of Entomological Research
http://dx.doi.org/10.1079/BER2003267
Service M.W., 1980, A guide to Medical Entomology, Macmillan International College Editions, London
Service M.W., 1985, Anopheles gambiae: Africa's principal malaria vector (1902-1984), Bulletin of the Entomological Society of America, 31: 8-12
Service M.W., 1988, New tools for old taxonomic problems in blood sucking insects in Biosystematics of haematophagous insects. ed. (M.W. Service), Clarendon press box, 325
White G.B. and Muniss, J.N., 1972, Taxonomic value of spermatheca size for distinguishing four members of the Anopheles gambiae complex in East Africa, Bull WHO, 46: 793-799
White G.B., 1977, The place of morphological studies in the investigation of Anopheles species complexes, Mosquito Systematics, 9: 1-24
WHO 1975, Manual on Practical Entomology in malaria: part II, Methods and techniques, WHO offset publications, No: 13. Geneva, Switzerland
Zahar A.R.M. and Hills and Davidson G., 1970, An attempt to group fresh water species of the Anopheles gambiae complex by some morphological larval and adult characters, Parasitologia (Rome), 12: 31-46